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The effect of density gradients on the supersonic wall modes (acoustic 
modes) of a 2-D confined compressible shear layer were investigated using 
linear analysis. Due to the inadequacies of the hyperbolic tangent profile, the 
boundary layer basic profiles were used. First a test case was taken with the 
same parameters as in Tam and Hu’s 15 analysis with convective Mach number 
M c = 1-836 and density ratio of 1.398. Three generalized inflection points 
were found giving rise to three modes. The first two show similar properties 
to the Class A and B modes in Ref. 15, and the third is an ‘inner mode’ which 
will be called a Class C mode. As the density ratio is increased, the smallest 
of the three neutral phase speeds tends towards the speed of the lower velocity 
stream, and the other two eventually coalesce and then disappear. These two 
effects lead to a linear resonance between the Class B modes which increases 
the cutoff frequency and growth rate of the lowest mode. In fact, growth rates 
of 2-4 times the test case were found as the density ratio was increased to 7. 
A similar trend is observed for the Class A modes when the density ratio is 
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decreased from the test case, but the growth rate is not changed by much from 
the test case. 
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I. Introduction 


The study of compressible shear layers has been thrust into the spotlight in 
the past decade mainly due to the need for supersonic combustors (SCRAM- 
JETS) in future generation aircraft. Numerous analyses have been conducted 
in order to better understand the instability characteristics and dynamics of 
compressible shear layers. In what follows a brief overview of some of these 
analyses will be presented with emphasis on work done using linear analysis. 
In the overview the term ‘constant density’ will refer to cases when the den- 
sity and temperature of the basic flow are taken to be constant throughout the 
shear layer so that their variation in the shear layer due to viscous heating is 
ignored. The researchers who have included viscous heating effects have used 
either the Busemann-Crocco (B-C) relations or some solution of the boundary 
layer equations for the density and temperature profiles. 

Some early work done on the inviscid instability of compressible shear lay- 
ers was reported in a series of papers by Blumen , Blumen et al. and Drazin 
and Garvey 3 . The work involved a 2-D temporal linear stability analysis of a 
compressible, constant density, unbounded shear layer modeled by a hyperbolic 
tangent profile based on the cross-stream coordinate, with equal and opposite 
velocities at the ends. For subsonic Mach numbers, they found a single mode 
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whose growth rate decreased as the Mach number increased. This mode is a 
Kelvin-Helmholtz (K-H) type of instability similar to that occurring in incom- 
pressible shear layers. As the Mach number is increased, the first order effect 
of compressibility is to decrease the growth rate of the 2-D K-H mode. For su- 
personic Mach numbers, two modes were found with equal but opposite phase 
speeds. These modes had very small growth rates compared to the subsonic 
K-H mode. Jackson and Grosch^ conducted a spatial stability analysis of a 
compressible shear layer (lower stream was quiescent) modeled by a hyperbolic 
tangent profile whose independent variable was a similarity variable based on 
Howarth’s transformation. Since this similarity variable is dependent on the 
Mach number, the hyperbolic tangent profile based on this transformation is 
also dependent on the Mach number. They computed the neutral modes for 
2-D and 3-D disturbances and also found two modes for supersonic Mach num- 
bers. One mode was supersonic with respect to the slow stream and subsonic 
with respect to the fast stream and was called a ‘fast mode’. The other was 
supersonic with respect to the fast stream and subsonic with respect to the 
slow stream and was called a ‘slow mode’. These modes were characterized 
by Mack 5 as radiating vorticity modes because of the oscillatory behavior of 
the eigenfunctions outside the chear layer. Sandham and Reynolds^ 6,7 ! inde- 
pendently confirmed the existence of these modes using a profile based on the 
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similarity solution of the compressible boundary layer equations. They also 

showed that for convective Mach numbers (Me), i.e. Mach number relative 

to the phase speed of the instability, above a critical subsonic value the most 

unstable disturbances were three-dimensional. This is in clear contrast to the 
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incompressible case. Similar results were obtained earlier by Gropengiesser 
and Lessen et Other studies on the unbounded compressible shear 

layer were conducted by Ragab and Wu^’^’^^, with the latter two references 
dealing with subharmonic instabilities in this type of shear layer. 

For later reference, it is also necessary to mention some results obtained 
by researchers concerning the use of the different profiles used for the above 
analyses. Sandham and Reynolds^ showed that use of the hyperbolic tangent 
profile based on the cross-steam coordinate was unsatisfactory for cases of 
hi gh shear or density ratio even at very low Mach numbers. They showed that 
errors of up to 16% in the growth rates can exist relative to the growth rates 
obtained by use of the boundary layer profiles. Jackson and Grosch 14 used 
three different profiles in investigating the fast and slow modes in a spatial 
shear layer. The first was a hyperbolic tangent profile based on the similarity 
variable mentioned above, the second was a profile based on the similarity 
solutions of the boundary layer profiles with a linear viscosity law and Prandtl 
number (Pr) equal to 1.0, and the third was a similar profile but with the 
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Sutherland viscosity law incorporated and Pr — 0.7. They showed that the 
results obtained by using the first profile were only qualitatively similar to the 
ones from the other two profiles. 

The analyses which more closely pertain to the current research were done 

on bounded (confined) shear layers as shown in Figure 1. Tam and Hu and 

Greenough et al conducted spatial and temporal analyses for a bounded 

shear layer. Greenough et al used a hyperbolic tangent profile based on the 

cross-stream coordinate and only considered the constant density case. They 

also found a single mode for subsonic Mach numbers. However, multiple modes 

were found when the Mach number was supersonic. These modes are different 

from the radiating vorticity modes found in the unbounded shear layer and 

were referred to as supersonic wall modes by Greenough or acoustic modes by 

Mack 5 . Tam and Hu 15 conducted a thorough spatial analysis of these modes 

using both a hyperbolic tangent profile (based on the cross-stream coordinate 

and using the B-C relations) as well as a vortex sheet model. For the vortex 

sheet model, they found two families of instability modes and labeled them 

Class A and Class B modes. They also found two families of neutral modes 

which they label Class C and D acoustic neutral modes. For the finite shear 

layer case, only the unstable solutions were shown because the contour of 
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integration was not deflected in their analysis to avoid the critical layer. 
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Tam and Hu found that for a relatively thick shear layer the most unstable 
mode is the first Class A 2-D mode. Mack^ discussed in detail the nature 
of these acoustic modes as they pertain to boundary layers, wakes, jets and 
confined shear layers. The underlying condition which needs to exist in order 
for acoustic modes to be observed is a region of trapped supersonic flow relative 
to the phase speed of the instability wave. This means that either two sonic 
lines (Me = 1.0), two walls or a combination of the two must exist in the flow. 
More results and facts about Tam and Hu’s analysis are given in the body of 
the present paper since they are needed for comparison to our own results. 

The effects of walls on the instabilities occurring in compressible shear 
layers were investigated by Zhuang et al and Morris et al . . They showed 

results indicating that the subsonic mode was not affected by the walls as 
long as the walls were not too close (10 vorticity thicknesses away). They also 
showed that, as the walls were brought in from infinity the radiating vorticity 
modes for supersonic Mach numbers first decreased in growth rate and then 
disappeared. Meanwhile, a new set of modes, which were not continuations of 
the vorticity modes, appeared and dominated the 2-D disturbances. 

There have been numerous numerical simulations done on unbounded shear 
layers such as Lele 20 , Sandham and Reynolds^ 6 ’ 7 ], Ragab and Wu’ 3 , and Ragab 
and Sheen 21 confirming the results obtained from linear analysis. Some have 
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also considered the bounded shear layer in their numerical simulations, such 
as Greenough et al ., Sigalla et Huang and Riley 2 **, Lu and Wu^ ’ ^ 

and Gathmann et a/. 27 . 

Many experiments have been conducted on the instability of shear layers, 
one of the most famous experiments being done by Brown and Roshko for 
very low subsonic Mach numbers. They showed the existence of large scale 
coherent structures which were formed within the shear layer and also showed 
that the growth of the shear layer depended linearly on the velocity ratio pa- 
rameter ( %+% )« They also looked at the effect of density ratio on the growth 
of the shear layer by taking three points = 7, 1, y) and observed that when 
the density of the lower(slower) stream is greater than the upper (faster) stream 
the growth rate increases from the uniform density case. The opposite was true 
when the upper stream was heavier than the lower one. Recently experiments 
on the compressible mixing layer have been conducted by Papamoschou and 
Roshko 29,30 . They observed a reduction in growth rate as a function of Mach 
number. They also showed that the growth rate curves could be parameterized 
by the use of a single Mach number which they defined as the Convective Mach 
Number ,M C . The definition of this Mach number is based upon using a refer- 
ence frame associated with the dominant wave. They defined two convective 


Mach numbers: 



( 1 . 1 ) 


M e i = 


Ui-Uc 
a l 


and 


M„ = d-2) 

a 2 

where U c is the speed of the dominant structure and ai,a2 are the speeds of 
sound in the two streams. By using an argument that the two streams share a 
common stagnation point, they showed that if the two streams have the same 
specific heat ratio then M*i = M c2 . Also from these arguments, it follows that: 


jj _ a 2^1 + a 1^2 

a\ + a2 


( 1 . 3 ) 


and 


M c = 


Ui^U2 

a l + a 2 


( 1 . 4 ) 


For small Mach numbers and almost equal specific heat ratios, they showed 
that one could express U c in terms of the density ratio: 


9 



( 1 . 5 ) 



Papamoschou'** has recently given a revised version of the above model 
for the case when the Mach number is high enough to produce shocks within 
the flow. In these cases it was found that M ci ^ M e 2, and so an asymmetric 
situation arises. This was predicted in the analysis of Jackson and Grosch in 
regard to their slow and fast modes, (i.e. the supersonic modes don’t obey 
the U c formula given above even for the constant density case). Papamoschou 
gave a revised model accounting for shock formation in the stream which is 
supersonic with respect to the dominant structure. More will be said regard- 
ing this result obtained by Papamoschou in later sections. The concept of 
convective Mach number has been applied to linear analysis by replacing Uc, 
the speed of the dominant structure, by c r , the phase speed of the instability. 

In this paper, first the governing equations and solution procedures will 
be discussed. The profiles of the basic flow velocity and density will be cal- 
culated from the boundary layer equations rather than assumed as in Tam 
and Hu 15 . The results for a test case matching the parameters in Tam and 
Hu’s analysis will then be presented. In the final portion, the results obtained 
when increasing and decreasing the density ratio from this test case will be 
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discussed. 


II. Flow Model and Governing Equations 

The model which is studied here is a confined compressible shear layer 
formed by two gases with different velocities, densities and properties, but 
with the same constant pressure. Figure 1 shows the configuration used in this 
analysis. The subscript 1 is used for the quantities related to the high speed 
freestream and the subscript 2 for the quantities of the low speed freestream. 
The stroamwise coordinate is x, the spanwise coordinate is z, and the cross- 
stream coordinate is y. For simplicity, free slip wall boundary conditions are 
assumed at the walls of the channel. It is also assumed that the flow is inviscid, 
non-conducting, and non-diffusive. For this situation the governing equations 
are the Euler equations for a two species system and in dimensional form can 
be written as the following: 


V 

dt* 


+ V-(pY) = 0 


( 2 . 1 ) 



+ v* • Vv*) + Vp* 


= 0 


( 2 . 2 ) 
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(2.3) 



4- v * • Vh*) 


dp* 


dt* 


- tf*. Vp* = 0 


d£g 

dt* 


+ tf-VC 0 =0 


(2.4) 


p* = p*R*T* (2.5) 

c 0 = 4. c*-4. p* =/>:+/>.* c„+c 6 =i (2.6) 

p* p* 

where p* is the density of the gas mixture, T* is the temperature, p* is the 
pressure, v* is the velocity vector, h* is the enthalpy which can be related to 
the temperature by h* = [c* Pa C a + c* pb ( 1 - Ca)]T* where c* Q and Cp b are the 
specific heat ratios at constant pressure, and C a is the mass fraction (species 
concentration) of species a. Also, R* = [i?aCa + i?j(l — Ca)\ is the gas constant 
of the mixture. 

These equations are non-dimensionalized using the fast (upper) freestream 
quantities, p * , U*, T* and the height of the channel, H*. Thermodynamic 
properties are also non-dimensionalized by the upper freestream thermody- 
namic properties, c* Pa and R* a . Based on this non-dimensionalization, the den- 
sity ratio ^ and the velocity ratio jj± are defined as p 2 and U 2 , respectively. 
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Once the equations are non-dimensionalized, they are linearized around a par- 
allel basic flow ( p(y ), P , U (y), Ca(y ))• These basic flow quantities can be found 
in several ways which will be discussed later. The flow variables can be written 
as: 


P- P(y) + p(x,y,z,t) (2.7a) 

p = P + p(x,y,z,t) (2-7 b) 

u = U(y) + u(x,y,z,t) (2.7c) 

v = v'(x,y,z,t) (2.7 d) 

w — w'(x,y,z,t) (2.7e) 

r = f(y) + T'(:r,y,M) (2-7/) 

Ca — C a {y) + c a (x, y, z, t) (2.7y) 
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Normal modes are assumed for these infinitesimal disturbances with the 


form: 


q = q(y ) exp[i(kx + flz — art)] (2-8) 

where q is the eigenfunction, k and (3 are the streamwise and spanwise wave 
numbers and u> is the frequency. In general fc, /3 and ui are complex. Once these 
normal modes are substituted into the equations of motion, a single O.D.E. 
can be found for the disturbance pressure eigenfunction: 



2kDU 
(w - kU) 


lDp}Dp+{^Y^-(u-kU) 2 -k' 2 -fi 2 }p = 0 (2.9) 


where 



- _ ^ 0 ^*° + ~ 

^ Cv a Ca + C»k(l — Ca) 

and Mi and 71 are, respectively, the Mach number and the ratio of specific 
heats of the high speed freestream. 
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The 2-D equation is obtained by setting = 0 (3-D disturbances will be 


considered later). 


D 2 p + { . 2kD ^,, - Id?}Dp + - ktlf - k 2 \p = 0 (2.10) 

t (u ) — kU) p > K 1 

The above equation is solved subject to the boundary conditions 

Dp = o y = ± \ ( 2 - n ) 

which come from the y-momentum equation by setting the normal velocity at 
the wall to zero. The rest of the eigenfunctions can be found from the following 
equations: 


Dp 

tp(u> — kU) 


(2.12 a) 


ikp + pDUv 
ip(u> — kU ) 


(2.126) 


P = 


p(ikU + Dv) + Dpi) 
i(C) - kU) 


(2.12c) 
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(2.12 d) 


vDCg 

i(u - kU) 


III. Basic Flow Profiles 

Before the method by which the eigenvalues axe obtained is discussed, it is 
necessary to discuss the basic velocity and density profiles used in this anal- 
ysis. The choice of profiles has varied from using piecewise constant profiles 
to profiles obtained from numerical simulations. One of the simpler ways 
adopted by many researchers is to assume that the basic velocity has a form 
of a hyperbolic tangent profile. As stated in the Introduction Sandham and 
Reynolds 6 have shown that this is not a good approximation for the actual 
velocity profile obtained by solving the compressible boundary layer equations 
when there is high shear or a high density ratio. A better approximation is the 
use of a scaled variable based on Howarth’s transformation for the variable in 
the hyperbolic tangent profile. Jackson and Grosch 14 demonstrated that this 
profile only gives qualitatively similar results when compared to the results ob- 
tained by using profiles found from the similarity solutions of the compressible 
boundary layer equations. Blanche and Reynolds showed that the profiles 
obtained by solving the O.D.E.’s obtained from the similarity transformation 
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of the compressible equations agreed very well with those found by numerical 
simulations. Because of the above reasons and reasons which will become ev- 
ident later in this section, the profiles chosen for this analysis were obtained 
by solving the compressible boundary layer equations after making use of a 
similarity variable. 

33 

The boundary layer equations for a multi-species flow are given in Kuo . 
For a two-species non-reacting system at constant pressure and Pr = Sc = 
Le = 1, they can be written as: 


dx 




(3.1a) 


♦ *du* , * *dv* d ( *du* \ 

p u 11= + p v or = W> 


(3.16) 


+ dh* f * +dhf d ( ♦ dhf \ 

“ dy* ~ dy* dy*> 


(3.1c) 


* *dCa * *dCa _ d / * dCg \ 
P U dx* P V dy* dy* dy* / 


(3.1d) 


* 

V 



(3.1e) 
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where h* = c*T* + - y- 

The solution procedure used for these equations is almost identical to that 
used in Planche and Reynolds with one difference being that in our analysis 
the product of the viscosity, p*, and the density, p*, is assumed to be constant 
(pp = 1) since the gas constant R is variable, instead of £ = 1 as in their paper. 
Also, the transformation from the similarity variable, r ;, to the cross-stream 
coordinate, y, is done by integrating the following equation: 


dy _ 1 I 2x 
drj p\ Re 


y(0) = 0 


and the appropriate value of the quantity y^r| is found from the condition 

Sfjj — 1. 

Finally the profiles are re-non-dimensionalized using the height of the chan- 
nel. In all cases considered the ratio of the vorticity thickness to the channel 
height is taken as ^ = 0.1 which corresponds to that used by Greenough et 
a/. 16 , Sigalla et a/J 22,23 ^ and Tam and Hu. This value of the vorticity thick- 
ness insures that the profiles have reached their freestream values far from the 
confining walls. It should be noted that the profiles used in this analysis are 
applicable for only small distances from the splitter plate. 

As stated in the Introduction, the test case that is taken here is the one 
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used in Tam and Hu’s analysis with the following properties: 


Mi = 4.5 M 2 = 1.6 


— = 1.29 7i = 1.67(Fe) 72 = 1A0(N 2 ) 

a 2 

Figure 2a, b and c show the profiles used in Tam and Hu’s analysis and the 
ones calculated from the boundary layer equations. Besides the quantitative 
differences seen in the velocity and density profiles, a major qualitative differ- 
ence is the existence of three extrema in the product of the basic density and 
the basic vorticity profiles pDU , as shown in Figure 2c. This phenomenon is 
dependent on the convective Mach number and the ratio of the specific heats. 
From Lees and Lin 34 , it is known that the zeros of this product can give rise 
to neutral modes. In the unbounded case, Jackson and Grosch 4 showed that 
only one of these extrema is related to a critical neutral mode. This is the 
subsonic neutral mode whose eigenfunctions decay exponentially away from 
the mixing layer. However, for the bounded case, all three extrema can be 
associated with critical neutral modes, as we shall show. 

It is also important to discuss the effect of the density ratio on the profiles 
which are calculated for our linear analysis. Figure 3 shows the velocity profiles 


19 



at a fixed convective Mach number and velocity ratio for different density 
ratios. It can be seen that as the density ratio is reduced from unity the high 
speed side has a “fuller” profile and the low speed side has a bigger defect. 
The opposite is true when the density ratio is increased from unity. Figure 4 
shows the effect of the density ratio on pDU. It is evident that the emergence 
of the three extrema at a fixed convective Mach number and velocity ratio 
is dependent on the density ratio. For the parameters considered here the 
three extrema exist for approximately 1.1 < P 2 < 1-7. The consequences 
of the multiple extrema will become evident in the later section where the 
characteristics of the instability modes are presented. 
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IV. Numerical Procedure For Finding The Eig envalues 


The numerical procedure used to find the eigenvalues is very similar to 
the method used in Tam and Hu. The method is a combination of the grid 
search method used for their vortex sheet problem and the numerical dispersion 
relation used for their finite shear problem and is described below. 

Outside the shear layer (|y| >> Du = 0. Therefore, equation (2.10) 
becomes the familiar vortex sheet equation with solutions: 

p(y) = A cos[Ai(l - y)] 6 <y <1 

p(y) — B cos[A 2 (l Ay)] — 6 > y > —1 

where 


A, = [M?(w - kV,f - fc 2 ] = 

(4.2a) 

= m 2 (4)% - *^) 2 - *¥ 

' a 2 ' 

(4.26) 


Following Tam and Hu, equation (2.10) is integrated by a fourth order 
Runge-Kutta method from -6 (below the shear layer where equation (4.1b) 
holds) to +6 (above the shear layer where equation (4.1a) holds). The nu- 
merical solution at y = 6 , say Bf(6), and its derivative are matched to the 


(4.1a) 

(4.16) 
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corresponding solutions of equation (4.1a), thus giving rise to the following 
eigenvalue relation for |6| >> jp. 

W — — cos[Ai(l — ^)]/ , (5) + Ai sin[Ai(l — 6)]f(6) — 0 (4.3) 

The above equation is zero for specific values of (o>, k). There are two view- 
points taken when solving this eigenvalue problem, namely, the spatial and 
temporal viewpoints. In a temporal viewpoint, the wavenumber, &, is taken 
to be real and the frequency, u>, is the complex eigenvalue which is found by 
solving the above eigenvalue problem. Conversely, in a spatial viewpoint, the 
frequency is taken to be real and the wavenumber is the complex eigenvalue. 
The viewpoint chosen for a particular system usually depends on whether the 
instability is absolute or convective in nature; a temporal viewpoint is adopted 
for absolute instabilities, and a spatial viewpoint is chosen for convective in- 
stabilities arising due to an upstream wavemaker. Rigorous mathematical 
techniques exist based on the behavior of the impulse response of a system 
which dictate the appropriate viewpoint. For detailed discussions on this sub- 
ject, the reader is referred to Briggs 35 , Bers 36 and Huerre and Monkewitz 37 . 
For the current analysis a spatial viewpoint is adopted since both streams of 
the shear layer are supersonic and therefore all the eigenvalues are initially 
taken to be located in the upper-half of the complex wavenumber plane where 
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the frequency is complex and its imaginary part is a large positive number 
(Tam and Hu 1 ^). 

To find the sets of values (u/ r , fcr+ifc* ) for which equation (4.3) holds, a grid 
search procedure is employed. A coarse grid is used on the domain of interest 
in the complex A:-plane. For a given u>, Im(W)=0 and Re(W)=0 are plotted 
by a contouring routine, and their intersection points are taken as the initial 
location of the roots. Then by the use of a Newton’s iteration routine these 
eigenvalues are further refined. When performing this procedure, one must 
keep two things in mind. First, when w and k are real, the integration contour 
for finding f(6) must be deformed into the complex y- plane to avoid the critical 
layer. Thus the contour of integration must be deformed in such a way that 
it is below the critical layer singularity but above the singularities of the basic 
profiles in the complex y-plane. Second, one needs to distinguish between 
amplifying and evanescent waves. A good explanation for this procedure is 
given in Briggs and Tam and Hu. By following Briggs, if one needs to find 
the amplifying waves for a given real frequency fi, then one first starts with 
a complex frequency + tflj, where f lj is a large enough positive number 
such that all the waves corresponding to the response of the system for x < 0 
are in the lower half of the complex wave number plane and all the waves for 
x > 0 are in the upper half plane. Then the imaginary part of the frequency 
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is incremented to zero, while keeping the real part constant. The roots that 
drop below the real k-axis which were originally in the upper half plane are 
the amplifying waves for x > 0. Similarly, the roots starting from the lower 
half plane and crossing real k-axis represent amplifying waves for the x < 0 
response. 

V. Parameters and their Variation 

In order to isolate the effect of basic flow density gradients from other 
effects such a- shear and compressibility, the following method is proposed. 
The convective Mach number based on the current non-dimensionalization is 
given by: 


Mr = 


U* - U* Afi(l - U 2 ) 


a? 


+ G* 


+ \f' 


- 32 - 

7iA>2 


(5.1) 


In order to remove the first order compressibility effects, it is chosen to 
work with a constant convective Mach number. The above convective Mach 
number was shown by Papamoschou^ to be a good parameter to describe the 
effects of compressibility. Also the velocity ratio = I7 2 is also taken to be 
constant. Jackson and Grosch 14 show that the maximum temporal growth 
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rate is proportional to (1 - U 3 ). Although the consequences of varying this 
term for spatial calculations has not been investigated it was decided not to 
vary this term for our purposes. Another way to think about this condition 
comes from recognizing that if this term is unchanged then the difference in 

m 

mean shear ( U * — U*) is unchanged. Therefore, as the density ratio = p 2 is 
changed, the fast stream Mach number, Mi, is changed to keep the convective 
Mach number constant. 

VI. Mode Labels 


The convention used r or the instability modes is the same as that used in 
Tam and Hu. All unstable modes whose phase velocity start from the velocity 
of the high speed stream will be labeled as Class Amn modes and all unstable 
modes whose phase velocity start from the velocity of the low speed stream 
will be labeled as Class B mn modes. The m and n subscripts are the 3-D 
and 2-D mode numbers respectively. The bulk of this paper will deal with the 
m = 0, i.e. a 2-D case. 

VII. Results and Discussion 
A. Test Case = 1.398 
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As stated earlier, the first case which is considered in this analysis is a 
test case based on the parameters used in the Tam and Hu analysis which 
correspond to Ale — 1.836, p? — 1.398, and U 2 = 0.276. For this case, the 
basic profiles have three generalized inflection points ( D(pDU ) = 0) which 
are summarized in Figure 5. Figure 6 shows the phase speed of the three 
different types of instability modes which are present for this case. The modes 
represented by the solid lines are the Class A modes which start from Cp/j 

~ JJ * which is a non-inflectional neutral solution. This type of solution 
k r 

is described in MadJ 5,38 ! in conjunction with the boundary layer acoustic 
instabilities and is possible only when a trapped region of supersonic flow 
relative to the phase speed of the instability wave exists in the flow. Some of 
these modes, such as A 01 , have wave velocities that are bounded by the velocity 
at the upper inflectional neutral point, and the others, such as Aq2 ? remain 
unstable. The modes represented by the dashed lines are modes which start 
with a wave velocity equal to the mean velocity at the middle inflectional 
neutral point and which apparently have not been discussed so far in the 
literature. These modes will be labeled as Class C modes. Some of these 
modes have wave velocities which are bounded by the velocity at the upper 
inflectional neutral point and the other remain unstable. The third set of 
modes which are represented by the dash-dot lines are the Class B modes. 
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These modes start from C p h = U 2 which is a non-inflectional neutral solution. 
The phase velocity of all of these modes are bounded by the velocity at the 
lower inflectional neutral point. It is clear that the use of a regular hyperbolic 
tangent profile based on the y coordinate would not give the results presented 
above (cf. Fig. 2c). Runs were made using a hyperbolic tangent profile, with 
the results being that both Class A and B modes terminated at the one and 
only inflection point. These results are shown in Figure 7. 

Figure 8a, b and c show the spatial growth rates of the Class A, B and C 
modes respectively. The interesting thing to note here is that the C 01 mode has 
the highest growth rate as opposed to the A n mode when the hyperbolic tan- 
gent approximation is used. Furthermore, these modes are associated with a 
local minimum of pDU . In constant density or incompressible flows, instability 
is, of course, associated only with a maximum of vorticity. 

The disturbance pressure eigenfunctions of the three modes at the maxi- 
mum growth rates are given in Figures 9a, b and c. If one defines two convective 
Mach numbers Me 1 and Mci as: 


M c 1 = 


U* - c* ph 


a* 


\/T — °P h 

Me 2 


dL - u: 


CL^ 


(7.1 a) 


(7.1 b) 
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where M c 1 is the convective Mach number for stream 1 and Afcj is that for 
stream 2 , then the modes with the maximum growth rates in each Class, i.e. 
Aoz,Bq4 and Cfoi, can be described in the following manner. The A03 mode 
is supersonic with respect to the slower stream and subsonic with respect to 
the high speed stream. This can be seen in the shape of the eigenfunction 
which is monotone in the upper half of the channel and oscillatory in the lower 
half. The opposite is true for the Bq^ mode. The C01 mode is supersonic 
with respect to both streams and shows an oscillatory behavior over the entire 
channel height. 

Another difference observed when using the profiles based on the compress- 
ible boundary layer equations as opposed to the hyperbolic tangent profiles is 
the magnitude of the eigenfunctions. When the hyperbolic tangent profiles are 
used the eigenfunction with the maximum magnitude is the streamwise pertur- 
bation velocity eigenfunction. This is not true for the boundary layer profiles 
where the perturbation density eigenfunction possesses the largest magnitude. 

B. Increasing Density Ratio 

Next, the density ratio is increased to p 2 = 3 . 0 . Different results were 
obtained for this case ranging from the change in basic profiles to the growth 
rates and eigenfunction magnitudes. First, the results for this density ratio 
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will be given and then the changes that occur between 1.398 < p 2 < 3.0 will 
be discussed by considering some intermediate cases. 

The basic profiles have only one generalized inflection point for p 2 = 3.0. 
The velocity at that point is close to that of the low speed stream. This is 
illustrated in Figure 10. Figure 11 shows the phase speed of the instability 
waves present in this case. The Class A modes show similar behavior to 
the ones present in the case of p 2 = 1.398; however there are more Class A 
modes present in the same frequency range considered. Also, since two of the 
inflectional neutral points are not present, these modes now remain unstable. 
The Class B modes, however, are very different. In this frequency range only 
the Bq\ mode has a growth rate which is higher than 0.01. The phase speeds 
of the rest of the Class B modes (not shown) start at the velocity at the non- 
inflectional neutral point Cph — U 2 and end at the velocity at the inflectional 
neutral point. These modes have very small growth rates. It will be shown 
later that, due to several factors as the density ratio is increased, a linear 
resonance occurs between the Class B modes which results in the particular 
configuration of the Class B modes shown in the figure. The growth rate of the 
Class A modes and the Bq i modes are given in Figures 12a and b. The many 
peaks present in the growth rate curve of the Bq\ mode is characteristic of the 
linear resonance (see below) which is occurring. It is interesting to note that 
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the peaks located at u> = 5.775 and u> — 6.62 have growth rates — kj = 0-7401 
and 0.7404 respectively. These two are almost identical and so both frequencies 
might be excited in the nonlinear portion of shear layer evolution. Also, this 
mode possesses the highest growth rate among all the modes and has twice 
the growth rate of the maximum growth rate for the p 2 = 1-398 case. This 
trend continues at higher density ratios and will be discussed further when the 
results for the case of p 2 — 7.0 are presented. 

Figure 13 shows the phase speed of the i?oi mode ^ can seen 

that the phase speed has an oscillatory behavior and becomes neutral after 
many oscillations. This behavior can also be attributed to the linear resonance 
and will also be discussed later. Another interesting result is observed in the 
behavior of the eigenfunctions. The eigenfunctions for the density and species 
concentration show large spikes near the critical layer and the magnitude of 
the density eigenfunction is two orders of magnitude higher than the velocity 
and pressure eigenfunctions. Figure 14a,b,c,d and e show the eigenfunction for 
the Boi mode at the first growth rate peak (w = 5.775) . The eigenfunctions 
at the second peak are very similar to those of the first peak. For comparison, 
Figure 14c also shows the perturbation density eigenfunction at the maximum 
growth of the test case which is normalized with respect to the eigenfunction 
magnitudes of the test case. 
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The next step will be to show the emergence of the different phenomena 
described above by analyzing the instability characteristics for density ratios 
between 1.398 and 3.0. As the density ratio is increased from 1.398, the upper 
and middle inflection points begin to move towards one another and the lower 
inflection point moves lower such that the velocity at that inflection point tends 
towards U 2 . At p 2 « 1.71 the upper and middle inflection points coalesce, and 
beyond that density ratio they disappear. This movement of the inflection 
points, i.e. changes in the basic profiles due to the effects of the density ratio, 
lead to a linear resonance between the Class B modes, as discussed below. 

In order to show that the interaction between the modes is indeed a res- 
onant interaction, a series of figures are provided showing how the modes 
coalesce and interact. First, Figure 15 shows the phase speed of the Class B 
modes for p 2 = 1.398 and p 2 = 1.8. As the density ratio is increased in this 
range the lower inflectional neutral phase speed decreased from 0.4 to 0.37. As 
this happens the curved parts of the phase speed curvse of the £oi mode and 
the B 02 mode connect at a density ratio between 1.398 and 1.8. As the density 
ratio is increased beyond that point, calculations show that the neutral phase 
speed decreases and the curves separate leaving the £ 0 1 the upper curved 
portion of the old £02 • This process continues and for p 2 = 1.8, as the figure 
shows, five of the original Class B modes have interacted. This type of process 
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is well described in Craik(1985). Figure 16 shows the growth rate curves for 
the two density ratios, and again it is possible to see the how the curve for 
pi — 1.8 evolves from the others. We should also note that for P 2 — 1*8 the 
growth rates of the Bq 2 ~ Bq§ modes are not shown because the magnitude of 
the growth rates is less than 0.01. Finally, Figure 17a and b show the phase 
speed and growth rates of the Class B modes for only pi — 1.8. On these 
figures the location where the higher modes become neutral are shown with 
the vertical dotted lines. The purpose of these figures is to show that these 
neutral frequencies correspond well with frequencies corresponding to the local 
minima of the Bq\ mode, cf. Fig. 16. 

Figure 18 shows the growth rate of the Bq\ for different density ratios. It 
is clear from this figure that as the density ratio is increased more and more 
Class B modes have resonant interactions with one another to cause the Bqi 
mode to have more and more peaks. Figure 19 shows the number of peaks in 
the growth rate curves as a function of the density ratio. It is clear that this 
is a non-linear step function of the density. 

Also, as the density ratio is increased, the phase speed of the Bq\ mode, 
which is the mode with the maximum growth rate, becomes more and more 
asymmetric. This type of phenomenon is observed in the experiments con- 
ducted by Papamoschou 31 . In his experiments, he observed that: “in su- 
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personic combinations M c i is always highest ( Uc —* U 2 ) while in supersonic- 
subsonic combinations, M e2 is always highest (Uc — 4 U\)” . In our analysis, an 
increase of the density ratio from the test case increased M 2 a-nd lowered M\ 
causing both Mach numbers to be between 2 and 4. Decreasing the density 
ratio increased M\ and decreased M 2 such that M 2 became almost sonic. 

The final comment which needs to be made concerns the effect of the 
density ratio on the growth rate of the modes. After all, this is one of the 
limiting factors in the design of supersonic combustors. Figure 20 shows the 
growth rate of the 5oi mode for p 2 = 7.0. This growth rate is nearly 400% 
greater than the growth rate at p 2 = 1.398. This type of increase in the growth 
rate as a function of the density ratio is supported by observations made by 
Brown and Roshko 28 who stated that growth rates in supersonic shear layers 
are more sensitive to changes in the density ratio than subsonic shear layers 
and that increases of up to 300% in growth rates were observed in supersonic 
shear layers. 

Figure 21 shows the maximum growth rate as a function of density ratio 
for the cases where mode £01 has become dominant. Density ratios from 2.0 
to 13 have been considered. Some leveling off of the curve is noticed as the 
density ratio is increased. 

The effect of the density ratio on the Class A modes varies from mode to 
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mode. For density ratios where Class C modes still exist, increasing the density 
ratio slightly increases or decreases the growth rate of some Class C modes and 
some Class A modes. These modes eventually coalesce and interact, and once 
the upper two inflection points disappear, only the Class A modes remain. As 
the density ratio is increased from 3.0 to 7.0 so that no Class C modes exist, 
the most unstable mode, Aoi, practically has a nearly constant growth rate, 
and the other Class A modes actually have lower growth rates. 

C. Decreasing the Density Ratio 

When the density ratio is decreased from the test case, almost the oppo- 
site of the previous case occurs. First, the lower and middle inflection points 
move towards one-another and eventually coalesce and disappear. The upper 
inflection point remains and tends upwards so that the inflectional neutral ve- 
locity tends towards U \ . These two effects should lead to a resonant interaction 
between the Class A modes. 

To investigate the above conjecture, two density ratios were considered, 
p 2 — 1 a nd p 2 = 7. In order to keep the convective Mach number constant, 
the gases in the mixing layer were changed to Ar in the fast stream and CH 4 
in the slow stream. This was done so that the Mach number of the slow stream 
remained supersonic (M 2 — 1.159 for P2 — 3 and M 2 — 1.000 for P2 = 7)* 
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Figure 22 shows the phase speed of the Class A and B modes for p 2 = 

The Class B modes ( Bq2 <md on) remain unstable, and within the same 
frequency range (0 < < 5) there are more Class B modes than in the 

p 2 = 1.398 case. A difference arises in the behavior of the Boi mode. For 
this density ratio, this mode becomes neutral at the upper inflectional neutral 
point. Recall that for p 2 = 3.0, the Aoi mode did not have the corresponding 
behavior. Now, the Aoi mode for this case interacts with the modes A03 — A05. 
However, the Aq2 mode is not involved in this interaction and remains unstable. 
Figures 23a and b show the growth rates for the Class A and B Modes. It is 
clear that both the A01 and Aq 2 modes have appreciable growth rates, and 
again like the higher Class B modes in the case of p 2 — 3.0, the higher Class 
A modes (A03 — A05) have growth rates which are less than 0.01. Due to the 
linear resonance occurring between some of the modes, the growth rate of the 
A01 mode (the maximum growth rate for this density ratio) is greater than 
the m aximum growth rate for the p 2 — 1.398 case. However, the maximum 
growth rate for this case is only slightly higher than the test case. The growth 
rate of the Class B modes are smaller for this density ratio as compared to the 
test case. 

When the density ratio is decreased to p 2 = 7, more interactions occur 
between the modes which is a bit more interesting. Figure 24 shows the phase 
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speed of the modes present at this density ratio. Comparing this figure to 
Figure 22 , one can see several changes. As the density ratio is decreased from 
^ to y, modes Bqi, Aqi and Aq 2 interact with one another which causes the 
I ?01 mode to be the dominant mode. Looking closely at the two figures the 
transition between them can be realized. First, just as there are more and more 
Class B modes in a certain frequency range as the density ratio decreases, there 
are fewer and fewer Class A modes. As the density ratio is decreased the Boi 
mode in Figure 22 interacts with the Aq 2 mode near the inflectional neutral 
point. Then this new Bqi mode, which has picked up the majority of the Aq 2 
mode’s curve, interacts with the A01 mode. This final interaction causes the 
Bq\ mode to be the mode which exhibits a linear resonance with the rest of the 
Class A modes. Actually at this density ratio more than 25 Class A modes are 
involved in this resonant interaction. Figure 25 shows the growth rate for the 
jBqi mode. Notice the unstable bandwidth which is many factors above any 
other mode considered at any density ratio. However, the increase in growth 
rate from the p? = 3 case is only about 0 . 05 . 

VIII. Concluding Remarks 

The effects of density gradients on the linear stability characteristics of 
the acoustic wall modes in a confined supersonic shear layer were studied. A 
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test case based on the analysis of Tam and Hu 15 was chosen. A new class 
of modes was discovered for this case which was due to the emergence of 
three generalized inflection points in the pDU profile. This new mode has the 
highest growth rate for the test case density ratio and is associated with the 
local minima of the pDU profile. 

Next, the case of increasing density ratio was considered. Increasing the 
density ratio caused the pDU profile to eventually have only one inflection 
point, with a velocity close to that of the low speed stream. This in turn 
led to a linear resonance between the Class B modes. As the density ratio 
was increased, more and more Class B modes were involved in the resonant 
interaction. Another phenomenon observed for this case was the increase in 
the growth rate as the density ratio increased. This increase in the growth 
rate seems to continue as the density ratio is increased with only some leveling 
of the curve noticed at the highest density ratio. Two effects can contribute to 
this increase. The first is the shape of the velocity profile. As the density ratio 
is increased the velocity is “fuller” on the low speed side which in turn results in 
higher growth rates as observed by Sandham and Reynolds in relation to the 
subsonic shear layer. The second effect has to due with the linear interaction 
of the modes which causes the “strengthening” of the first mode. 

As the density ratio was decreased, the pDU profile also lost two out of 
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the three inflection points. However, the remaining inflection point was the 
opposite of the one for the increasing density ratio case, with the velocity at 
the remaining inflection point being close to that of the high speed stream. 
This led to an interaction between the Class A modes, and for very low density 
ratios, also the Hoi mode. The increases in the growth rates for this case are 
small compared to the increasing density ratio. This is due to the fact that as 
the density ratio is decreased the velocity profile is “fuller” on the high speed 
side which tends to reduce the growth rate. This nearly offsets the effect of 
the linear resonance. 

A major conclusion which can be drawn from this analysis is that when 
a much heavier gas occupies the slow stream of a supersonic confined mixing 
layer, the growth rate can be increased by 4-5 factors from the case when 
the densities of the two gases are nearly equal. These high growth rates are 
comparable to the incompressible growth rates and could lead to better mixing 
in different applications. 


38 



ACKNOWLEDGMENTS 


This work was supported by NASA Grant NCC2-374. 


39 



References 


*Blumen, W. “Shear layer instability of an inviscid compressible fluid”, J. 
Fluid Mech., 40, 769-781, (1970). 

2 Blumen, W., Drazin, RG. and Billings, D.F. “Shear layer instability of 
an inviscid compressible fluid. Part 2”, J. Fluid Mech., 71, 305-316, (1975). 

^Drazin, P.G. and Garvey, A. “Shear layer instability of an inviscid com- 
pressible fluid. Part 3”, J. Fluid Mech. 82, 255-260, (1977). 

4 Jackson, T.L. and Grosch, C.E. (1989) “Inviscid spatial stability of a 
compressible mixing layer”, J. Fluid Mech., 208, 609-637, (1989). 

8 Mack, L. M., “On the inviscid acoustic-mode instability of supersonic 
shear flows-Part 1: Two-Dimensional waves”, Theoret. Comput. Fluid 
Dynamics, 2, 97-123, (1990). 

6 Sandham, N.D. and Reynolds, W.C., “A numerical investigation of the 
compressible mixing layer”, Report No. TF-45. Thermosciences Division, 
Mechanical Engineering Department, Stanford University, (1989). 

7 Sandham N.D. and Reynolds, W.C. “Three-Dimensional simulations of 
large eddies in the compressible mixing layer”, J. Fluid Mech., 224, 133- 
158, (1991). 

8 Gropengiesser, H. “Study on the stability of boundary layers and com- 
pressible fluids”, NASA TT F12786, (1970). 

9 Lessen, M., Fox, J.A. and Zien, H.M. “On the inviscid stability of the 
laminar mixing of two parallel streams of a compressible fluid”, J. Fluid 
Mech., 23, 355-367, (1965). 

10 Lessen, M., Fox, J.A. and Zien, H.M. “Stability of the laminar mixing 
of two parallel streams with respect to supersonic disturbances”, J. Fluid 
Mech., 25, 737-742, (1966). 

^Ragab, S.A. and Wu, J.L. “Linear instabilities in two- dimensional Com- 
pressible mixing layers”, Phys. Fluids A, 1, 957-966, (1989). 

12 Ragab, S.A. and Wu, J.L. “Linear subharmonic instabilities of periodic 
compressible mixing layers”, AIAA Paper No. 89-0039, (1989). 



13 Ragab, S.A. and Wu, J.L. “Instabilities of supersonic shear flow”, AIAA 
Paper No. 90-0712, (1990). 

14 Jackson, T.L. and Grosch, C.E. “Inviscid spatial stability of a compress- 
ible mixing layer. Part 3. Effect of thermodynamics”, J. Fluid Mech., 224, 
159-175, (1991). 

15 Tam, C.K.W. and Hu, F.Q. “The instability and acoustic wave modes 
of supersonic mixing layers inside a rectangular channel”, J. Fluid Mech., 
203, 51-76, (1989). 

16 Greenough, J., Riley, J., Soetrisno, M. and Eberhardt, D. “The effect of 
walls on a compressible mixing layer”, AIAA Paper No. 89-0372, (1989). 

17 Hu, F.Q. Private Communication (1994). 

18 Zhuang, M., Dimotakis, P.E., and Kubota, T., “The effect of walls on 
a spatially growing supersonic shear layer”, Phys. Fluids A, 2 599-604, 
(1990). 

^Morris, P.J., and Giridharan, M.G., “The effect of walls on the instability 
waves in supersonic shear layers”, Phys. Fluids A, 3 356-358, (1991). 

20 Lele, S.K. “Direct numerical simulation of compressible free shear flows”, 

AIAA Paper No. 89-0374, (1989). 

21 Ragab, S.A. and Sheen, S., “Numerical simulation of a compressible mix- 
ing layer”, AIAA Paper No. 90-1669, (1990). 

22 Sigalla, L., Eberhardt, D., Greenough, J., Riley, J. and Soetrisno, M. 
“Numerical simulation of confined, spatially-developing mixing layers: com- 
parison to the temporal shear layer”, AIAA Paper No. 90-1462, (1990). 

23 Sigalla, L., Eberhardt, D., Greenough J. and Riley, J. “Numerical sim- 
ulation of confined, spatially developing shear layers: Nonlinear interaction 
between instability modes”, AIAA Paper No. 91-1643, (1991). 

24 Huang, H.-S. and Riley, J.J. “Transition to turbulence in confined, com- 
pressible mixing layers. Part I: 3-D numerical simulations with excitation of 
random, broadband white noise”, AIAA Paper No. 92-05532, (1992). 

25 Lu, P.J. and Wu, K.C. “On the shock enhancement of confined supersonic 
mixing layers”, Phys. Fluids A, 3, 3046-3062, (1991a). 



26 Lu, P.J. and Wu, K.C. “Numerical investigation on the structure of a 
confined supersonic mixing layer”, Phys. Fluids A, 3, 3063-3079, (1991b). 

27 Gathmann, R.J., Si-Ameur, M. and Mathey, F. “Numerical simulations 
of three-dimensional natural transition in the compressible confined shear 
layer”, Phys. Fluids A, 5, 2946-2968, (1993). 

23 Brown, G.L. anc Roshko, A. “On density effects and large structure in 
turbulent mixing layers”, J. Fluid Mech., 64, 775-816, (1974). 

29 Papamoschou, D. and Roshko, A. “Observations of supersonic free shear 
layer”, AIAA Paper No. 86-0126, (1986). 

30 Papamoschou, D. and Roshko, A. “The compressible turbulent shear 
layer; an experimental study”, J. Fluid Mech., 197, 453-477), (1988). 

3 *Papamoschou, D. “A structure of the compressible turbulent shear layer” , 

AIAA Paper No. 89-0126, (1989). 

32 Plancht, O.H. and Reynolds, W.C. “A Numerical investigation of the 
compressible reacting mixing layer”, Report No. TF-56. Thermosciences Di- 
vision, Mechanical Engineering Department, Stanford University, (1992). 

33 Kuo, K.K., Principles of Combustion, Wiley Interscience (1982). 

34 Lees, L., and Lin, C.C., “Investigation of the stability of the laminar 
boundary layer in a compressible fluid”, NACA Technical Note No. 1115, 
(1946). 

35 Briggs, R.J. Electron-Stream interaction with plasmas, Research 
Monograph No. 29. Cambridge, Mass: M.I.T. Press, (1964). 

36 Bers, A. “Linear waves and instabilities”, In Physique des Plasmas 
(ed. C. DeWitt & J. Peyraud), pp. 117-213, (1975). New York: Gordon and 
Breach. 

37 Huerre, P. and Monkewitz, P.A. (1990) “Local and global instabilities in 
spatially developing flows”, Annu. Rev. Fluid Mech., 22, 473-537, (1990). 

33 Mack, L.M., “Boundary-Layer stability theory”. In Special Course on 
Stability and Transition of Laminar Flow, AGARD Report No. 709, 3-1 
to 3-81, (1984). 



List of Figures 


Figrue 1. Model used for the analysis of the compressible shear layer. 

Figure 2a. Base velocity profiles for the test case. 

Figure 2b. Base density profiles for the test case. 

Figure 2c. Base density times base vorticity ( pDU ) profiles for the test 
case. 

Figure 3. Base velocity profiles calculated from the boundary layer equa- 
tions for different density ratios at a fixed convective Mach number ( M c = 
1.836) and velocity ratio (17 2 = 0.276). 

Figure 4. Base density times base vorticity ( pDU ) calculated from the 
boundary layer equations for different density ratios at a fixed convective Mach 
number (M c = 1.836) and velocity ratio (U 2 = 0.276). 

Figure 5. Graph o* D(pDU) for the test case. Zeros are given by the 
vertical dotted lines. (Boundary layer profiles) 

Figure 6. Phase speed vs. frequency for the unstable modes present in the 
flow at the test case parameters. 

Figrue 7. Phase speed vs. frequency for the unstable modes present in the 
flow at the test case parameters. (Hyperbolic Tangent profiles) 

Figure 8a. Spatial growth rate vs. frequency for the Class A modes at the 
test case parameters. 

Figure 8b. Spatial growth rate vs. frequency for the Class B modes at the 
test case parameters. 

Figure 8c. Spatial growth rate vs. frequency for the Class C modes at the 
test case parameters. 

Figure 9a. Disturbance pressure eigenfunction for the Aq 3 mode at the 
maximum growth rate. (Test case parameters) 

Figure 9b. Disturbance pressure eigenfunction for the Bq 4 mode at the 
maximum growth rate. (Test case parameters) 



Figure 9c. Disturbance pressure eigenfunction for the Coi mode at the 
maximum growth rate. (Test case parameters) 

Figure 10. Graph of D(pDU ) for p 2 = 3.0. Only one zero exists and is 
given by the vertical dotted line. 

Figure 11. Phase speed vs. frequency for the unstable modes present in 
the flow at p 2 = 3.0. Higher Class B modes ( i?02 and on ) not shown. 

Figure 12a. Spatial growth rate vs. frequency for the Class A modes at 
p 2 = 3.0. 

Figure 12b. Spatial growth rate vs. frequency for the Hqi mode at p 2 = 3.0. 

Figure 13. Phase speed vs. frequency for the Hoi mode at p 2 = 3.0. 

Figure 14a. Streamwise disturbance velocity eigenfunction for the Hqi 
mode for p 2 = 3.0 and at the first maximum growth rate peak u > = 5.775). 

Figure 14b. Cross-stream disturbance velocity eigenfunction for the Hoi 
mode for p 2 = 3.0 and at the first maximum growth rate peak u> = 5.775). 

Figure 14c. Disturbance density eigenfunction for the Hoi mode for p 2 - 
3.0 and at the first maximum growth rate peak u> = 5.775), and for the Coi 
for p 2 = 1.398 at c o — 1.75). 

Figure 14d. Disturbance pressure eigenfunction for the Hoi mode for p 2 — 
3.0 and at the first maximum growth rate peak u> = 5.775). 

Figure 14e. Disturbance species concentration eigenfunction for the Hoi 
mode for p 2 = 3.0 and at the first maximum growth rate peak u> = 5.775). 

Figure 15. Phase speed vs. frequency for the Class B modes at density 
ratios of 1.398 and 1.8. 

Figure 16. Spatial growth rate vs. frequency for the Class B modes at 
p 2 = 1.398 and the Hoi mode at p 2 = 1.8. 

Figure 17a. Phase speed vs. frequency for the Class B modes at p 2 = 
1.8. Vertical dotted lines represent locations where the higher modes become 
neutral. 



Figure 17b. Spatial growth rate vs. frequency for the Bqi mode at p 2 - 
1.8. Vertical dotted lines represent locations where the higher modes become 
neutral. 

Figure 18. Spatial growth rate vs. frequency for the Boi mode at different 
density ratios. 

Figure 19. Growth rate peaks (number of oscillatior :>) vs. density ratio for 
the Bqi mode. 

Figure 20. Spatial growth rate vs. frequency for the Bq\ mode at p 2 = 7.0. 

Figure 21. Maximum spatial growth rate vs. density ratio for the Bq\ 
mode. 

Figure 22. Phase speed vs. frequency for the unstable modes present in 
the flow at p 2 = 3. 

Figure 23a. Spatial Growth rate vs. frequency for the Class A modes at 
„ - 1 
p 2 - 3. 

Figure 23b. Spatial Growth rate vs. frequency for the Class B modes at 

. - 1 
p 2 - 3. 

Figure 24. Phase speed vs. frequency for the unstable modes present in 
the flow at p 2 = y. 

Figure 25. Spatial Growth rate vs. frequency for the Bq\ mode at p 2 = j. 





OSO- 




.5 (He) U 2 = 0.276 






















Stream 1 He 







CO 

t^- 

CM 

• 

O 

It 

l\l 

ZD 


GO O 
O) O 

« ^ 

T— O 


II II 

Ovj 8 

CL tO 


O O CO 

o o CO 

ID CD 00 
• • • 

^ ▼“ 
II II II 


CM 



O 










• 111 ! 















.500 (He) U, 
.600 (N,) p a 
.836 5. 




















■ 4.500 (He) U 










4.500 




Figure 10 


























3.876 (Ho) 
2.019 (N 2 ) 
1.836 












6 (He) 








Figure 14b 










3.876 

2.019 



- 0.2 











= 3.876 (He) Uj = 0.276 






3.876 (He) 
2.019 (N,) 



eo- 








0.45 



0.00 0.75 1.50 2.25 







Figure 16 








2.25 3.00 




Figure 17b 










Figure 19 



e/uno 0jeu mimojo ui «M*®d 1° J®9 u,n N 


2.25 2.50 2.75 3.00 3.25 



Figure 20 



25.0 




Figure 21 




























